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Abstract 

Intermediate energy (p, p'x) reaction is studied with antisymmetrized molec- 
ular dynamics (AMD) in the cases of 58 Ni target with E p = 120 MeV and 
12 C target with E p = 200 and 90 MeV. Angular distributions for various E p * 
energies are shown to be reproduced well without any adjustable parameter, 
which shows the reliability and usefulness of AMD in describing light-ion re- 
actions. Detailed analyses of the calculations are made in the case of 58 Ni 
target and following results are obtained: Two-step contributions are found 
to be dominant in some large angle region and to be indispensable for the re- 
production of data. Furthermore the reproduction of data in the large angle 
region > 120° for E p > = 100 MeV is shown to be due to three-step contri- 
butions. Angular distributions for E p > •> 40 MeV are found to be insensitive 
to the choice of different in-medium nucleon-nucleon cross sections o~nn an d 
the reason of this insensitivity is discussed in detail. On the other hand, the 
total reaction cross section and the cross section of evaporated protons are 
found to be sensitive to o~nn- In the course of the analyses of the calculations, 
comparison is made with the distorted wave approach. 
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I. INTRODUCTION 



Nucleon inelastic scattering to continuum at intermediate energies has been extensively 
studied, through which our understanding of pre-equilibrium processes in the scattering has 
much advanced. Among many kinds of theoretical investigations including the exciton model 
|l|fj] and other multi-step reaction theories , the intra- nuclear cascade model (INC) ||[7j 
has served as an important approach since the activation of many degrees of freedom in pre- 
equilibrium processes is complicated compared to the compound-nucleus reaction process 
with full equilibrization. Because the original INC model has several drawbacks such as the 
absence of the attractive interaction among nucleons, many kinds of modifications have been 
introduced into the INC model. 

Recently, besides INC and its modified versions, new types of transport theories (or mi- 
croscopic simulation theories) like BUU (Boltzmann-Uehling-Uhlenbeck) ||, QMD (quan- 
tum molecular dynamics) |5],|TU], and AMD (antisymmetrized molecular dynamics) JTT|JT^ | 
have been developed in order to investigate complicated reaction processes of heavy-ion col- 
lisions at intermediate and high energies. These new types of transport theories can describe 
the self-consistent mean field of the system which changes with time depending on the stage 
of the reaction process. These theories are, of course, also applicable to light-ion reactions 
including the nucleon inelastic scattering to continuum. 

We think that QMD and AMD approaches to reactions induced by light ions and also 
hadrons and leptons are especially important. These molecular dynamics models can de- 
scribe dynamical production processes of fragments and hence they can provide us with a 
unified theoretical description of two different kinds of reaction processes, namely processes 
of the type of light-ion physics and those of the type of heavy-ion physics. Nucleon inelastic 
scattering to continuum belongs to the former type and fragmentation reaction to the latter 
type. 

In this paper we study proton inelastic scattering to continuum at intermediate energies 
by the use of the AMD model. The AMD model is a new transport theory and has already 
proved to be very successful in describing heavy- ion collisions at medium energies fl2|-18 



AMD describes the total system with a Slater determinant of nucleon wave packets and 
hence it has quantum mechanical character, which has been demonstrated in the ability 
of treating shell effects in the dynamical formation of fragments. Furthermore it has been 
shown that ground-state wave functions of colliding nuclei given by the AMD model are 
realistic and reproduce many spectroscopic data very well |19| -2T|]. We report in this paper 
the study of (p, p'x) reaction in the cases of 58 Ni target with E p = 120 MeV and 12 C target 
with E p = 200 and 90 MeV. We show that the angular distributions are well reproduced 
by AMD without any adjustable parameter. This shows the reliability and usefulness of 
AMD in treating light-ion reactions. We make detailed analyses of the calculations in the 
case of 58 Ni target in the following way: We decompose the calculated cross sections into 
the contributions coming from different steps in order to study magnitudes of multi-step 
contributions. Two-step contributions are shown to be dominant in some large angle region 
and to be indispensable for the reproduction of data. Furthermore the reproduction of data 
in the large angle region 9 > 120° for E p r = 100 MeV is shown to be due to three-step 
contributions. Angular distributions for E p i 40 MeV are shown to be insensitive to the 
choice of different in-medium nucleon-nucleon cross sections u^n and the reason of this 
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insensitivity is discussed in detail. On the other hand, the total reaction cross section and 
cross section of evaporated protons are shown to be sensitive to <tnn- In discussing the 
calculated results we make comparison with the results obtained with the semi-classical 



distorted- wave approach of Ref. [p2 



The organization of this paper is as follows: In Sec. II we explain the AMD framework, 
the adopted effective two-nucleon force, and three choices of in-medium nucleon-nucleon 
cross section <Jnn- A detailed explanation of the definition of the step number of reaction 
process is also given. In Sec. Ill we give the comparison with experiments of the calculated 
angular distributions at various E p t energies. Here the decomposition of the cross sections 
into multi-step contributions is also made. In Sec. IV we make detailed analysis of the 
dependence of the calculation on the different choice of <tnn- Finally in Sec. V we give the 
summary. 



II. FORMULATION 



A. AMD 



The framework of AMD (antisymmetrized molecular dynamics) was described in detail 



in Ref. |12| and hence we here explain only the outline of the AMD theory. 

In AMD, the wave function of A-nucleon system is described by a Slater determinant 
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where \ stands for the spin-isospin function and ay represents the spin-isospin label of the 
j-th single particle state, aj — p t, p I, n |, or n J,, 0z. is the spatial wave function of the 
j-th single-particle state which is a Gaussian wave packet, 
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where the width parameter v is treated as time-independent in the present work. We take 
^=0.16 fm -2 in the calculation in this paper. Here Zj is the complex vector whose real 
and imaginary parts, Dj and Kj, are the spatial and momentum centers of the packet, 
respectively. 

The time developments of the coordinate parameters, Z = {Zj ■, (J = 1, 2, . . . , A)}, are de- 
termined by two processes. One is the time development determined by the time- dependent 
variation-al principle; 
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which leads to the equation of motion for Z, 
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A a ($(z)\Hmz)) 



(4) 



where a,r = x, y, z. 

During the dynamical reaction stage, the total system can be separated into several 
isolated nucleons and fragments. Since the wave functions of the center-of-mass motion of 
these isolated nucleons and fragments are Gaussian wave packets, each of these isolated 
particles carries spurious zero-point energy of its center-of-mass motion. The total amount 



of the spurious energy of center-of-mass motion can be expressed as a function of Z ||12| . |11 
which we denote as E sprs (Z). The actual Hamiltonian we use in the above equation of motion 
(Eq. (D) is, therefore, given by {Q(Z)\H\${Z)) / (${Z)\$(Z)) - E sprs (Z). 

The second process which determines the time development of the system is the stochastic 
two-nucleon collision process. We incorporate this process in a similar way as it is done in 
QMD by introducing the physical nucleon coordinates {Wj}, mimicking the time- dependent 
cluster model (TDCM) @, as 



i A 



(5) 



dm ■ z 



Here it should be noted that, due to the antisymmetrization, Dj and K,, of Zj do not 
represent the position and momentum of j-th nucleon, respectively. When physical nu- 
cleon positions Rj and R& become close each other, these j-th and fc-th nucleons can 
make stochasic two-nucleon collisions. Let initial Wj and be changed into final Wj 
and W h by a two-nucleon collision. In order to continue the calculation of time de- 
velopment of the system wave function after this collision, we need to back-transform 
{W 1; . . . , Wj, . . . , W k , . . . , W^} into {Z[, Z' 2 , . . . , Z' A }. However, in general, the back- 
transformation from W = {Wj (j = 1,...,A)} to Z = {Zj (J = I,..., A)} does not 
always exist. When the back-transformation does not exist, we regard that the two-nucleon 
collision is Pauli-blocked. W is defined to be in Pauli-forbidden region if it cannot be back- 
transformed to any Z. The notion of the Pauli-forbidden region defined above is an extension 
of that of TDCM ||23|| . AMD without stochastic two-nucleon collisions is the same as FMD 
(Fermionic Molecular Dynamics) |23]], and FMD is a special case of TDCM [p5| , p3[ where 
every cluster is composed of a single nucleon. 

The full procedure of the AMD description of the nuclear reaction consists of three major 
steps: The first step is the initialization, namely the construction of the wave functions of 
the ground states of colliding nuclei. The initialization is made by the use of the frictional 
cooling method [p~T| ? pT9| |2T|1 . It has been checked that wave functions given by AMD are 



realistic and reproduce many spectroscopic data very well. The second step is the calculation 
of dynamical collision stage by the equation of motion and stochastic two-nucleon collisions. 
The final step is the calculation of the statistical decay of primordial fragments. Primordial 
fragments mean the fragments which are present when the dynamical stage of the reaction 
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has finished. These fragments are not in their ground states but are excited, and they decay 
through evaporation with a long time scale. In this paper, the switching time from the 
dynamical stage to the evaporation stage was chosen to be 150 fm/c. Statistical decays of 
fragments were calculated with the code of Ref. [26| which is similar to the code of Puhlhofer 
271. 



B. Effective force and in-medium nucleon-nucleon cross sections 



As the effective two-nucleon force, we adopt the Gogny force [28|] which has been sue 



cessfully used in studying heavy-ion reactions with AMD []T4| , [T5|j . The Gogny force consists 



of finite-range two-body force and density-dependent zero-range repulsive force. This force 
gives a momentum-dependent mean field which reproduces well the observed energy depen- 
dence of the nucleon optical potential up to about 200 MeV. The nuclear matter EOS given 
by this force is soft with the incompressibility K = 228 MeV. Corresponding to the choice of 
Gogny force, the calculational formula of the total spurious center-of-mass energy E sprs (Z) 
is taken to be the same as Ref. |I4J in the case of 12 C target, while in the case of 58 Ni target 



the value of the To parameter in the formula of E sprs (Z) is changed into 8.70 MeV leaving 
other parameters unchanged. The binding energies of 12 C and 58 Ni are calculated to be 92.6 
MeV and 507.6 MeV, respectively, while the observed values are 92.2 MeV and 506.5 MeV, 
respectively. The r.m.s. radii of 12 C and 58 Ni are calculated to be 2.55 fm and 3.85 fm, 
respectively, which are reasonable. 

As the in-medium nucleon-nucleon cross section <Jnn, we adopt the following three dif- 
ferent ones, case-1, case-2, and case-3: Case-1 a^^ including its angular distribution is the 
same as the a^jy used in Ref. |1J and is given asffjvw = min {cr^y, p^n}- Here a§ N defines 



the cross section for the high energy region and is the same as the free cross section which 
is parametrized as 



a£ n = max{l3335(E/MeV)- L125 ,40} mb, 

= a* n = max{4445(£/MeV)- L125 ,25}mb. 



(6) 



where E is the laboratory energy. On the other hand, aj^ N defines the cross section for the 
low energy region. It is density-dependent and is given as 



L _ L _ L 
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100 mb 



1 + E/(200MeV) + Cmin {(p/p ) 1/2 , 1} ' 
C = 2. 



(7) 



where po is the normal density, po = 0.17 fm 3 . The center-of-mass angular distribution of 
pp and nn collisions is taken to be isotropic while that of pn collisions is given as 
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Case-2 a^N is the same as case-1 a nn except that C of a NN ls taken to be C = instead 
of C = 2. Hence case-2 <r NN is not dependent on density. Case-3 ctjvjv including its angular 
distribution is the same as the cross section adopted in Ref. [p2] . It is just the free cross 
section and its parametrization is taken from Ref. as follows; 

/34.10 82.2 \ , 

o vn = (-= r + 82.2) mb, 



P 2 P 

10.63 29.92 

T 2 J 



cr„ n = 7T- + 49.9) mb, (9) 



„ /2B 



c V Tnc 2 



The center-of-mass angular distribution is taken from Ref. [29], where pp and nn collisions 



is isotropic while that of pn collisions for 40 MeV < E < 280 MeV is given as 

dapn j A x (E) + B x (E) (cos 6) 3 , for < 9 < tt/2 
dfi ^ \ A X {E) + B 2 (E)(cos9) 4 , for tt/2 < # < vr [ Uj 

with values of B\(E), and B2(E) being shown in Table I. The pn collisions below 40 

MeV is taken to be isotropic. 

In this paper we study 12 C(p, p'x) by using only case-1 <tnn, while we apply all the three 
kinds of a NN to the study of 58 Ni(p, p'x). 



C. Step number of reaction process 

Here for the sake of later discussion of multi-step processes, we define the number of 
steps in our AMD approach. What we need is to determine for each out-coming proton the 
step number of the reaction process. The step number should be defined as the number of 
collisions which have contributed in emitting the nucleon. If an outcoming proton originates 
from the decay of some primordial fragment, the process which this proton has experienced 
is a compound-nucleus process and we need not to define the step number; namely the step 
number is concerned only with protons which are emitted dynamically. 

The step numbers are defined and calculated in the following way. We put a label to 
every nucleon so that the label of the incident proton is 1. Let us consider the first collision 
of the incident proton with a target nucleon with label k. After the first collision, each of 
the two nucleons 1 and k may further make collisions with other target nucleons. We put 
ordering numbers to all these successive collisions beginning with the first collision. We 
call the ordering number of collision the collision index. The collision index of the first 
collision is No.l. Just after the first collision, we let each nucleon i of the total system 
carry its respective set Cj(l) composed of related previous collision indices. The sets Ci(l) 
and Cfc(l) of the two nucleons 1 and k are Ci(l) = Cfc(l) = {No.l} but the sets Cj(l) of 
other nucleons ( % ^ 1, k ) are all null, namely Cj(l) = (z ^ 1, k). One of the nucleons 1 
and k can make the collision with collision index No. 2. Let us consider the case that the 
nucleon k makes the No. 2 collision with a nucleon j. Then, just after the No. 2 collision we 
let the nucleons k and j carry not the old sets Cfc(l) = {No.l} and Cj(l) = but new sets 
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Cfc(2) = Cj(2) = {No.l, No. 2}. However, the sets Cj(2) of other nucleons than k and j after 
the No. 2 collision are made unchanged, namely Cj(2) = Cj(l) for i 7^ A;, j. In general, if a 
nucleon m makes the No.p collision with a nucleon n, the new sets after the No.p collision 
are 

C m (p) = C n ( P ) = C m (p - 1) U C n (p - 1) U {No.p}, 
Ci(p) = Ci(p — 1) for i m,n. 

Note that the double counting of the same collision index is to be avoided in constructing 
new sets {Cj(p)} from old sets {Cj(p — 1)}. 

What is important in the above rule is that we only consider two-nucleon collisions which 
occur successively starting with the first collision of the incident proton as we explained above 
in time order. We give collision indices only to these two-nucleon collisions and we call them 
indexed collisions. Any other two-nucleon collisions which are not induced by the first 
collision of the incident proton have no collision index and are called non-indexed collisions. 
Non-indexed collisions play no role in constructing the sets {Ci(p)}. For example, the set 
Ci(p) of a nucleon i ^ 1 remains to be a null set for any p if this nucleon i is not involved 
at all in any indexed collision, even when it experiences many non-indexed collisions. 

Thus the set C$(p) is a set composed of all the indexed collisions that have had influence 
on the z-th nucleon at the moment just after the No.p collision. If the total number of 
collision indices contained in C{(p) is N, the nucleon % after the No.p collision is regarded as 
having a step number N. Especially when an outcoming nucleon has a step number N, this 
nucleon is defined to be due to iV-step process. 

In Fig. 1 we give two illustrative examples of multi-step process. In both cases (i) and 
(ii) of this figure, the nucleon (a) is the incident proton and collision points are labeled by 
collision indices 1, 2, . . . instead of No.l, No. 2, .... In the case of (i), both of outcoming 
nucleons (b) and (c) carry the same set of collision indices {1,2,3} and hence are three- 
step nucleons. In the case of (ii), the nucleon (b) carries a set of collision indices {1,2,3}, 
and hence is a three-step nucleon, while the nucleon (c) carrying a set of collision indices 
{1,2,3,4} is a four-step nucleon. We should note here that the three-step nucleons (b) and 
(c) in the case of (i) are formed as the result of a different type of collision process from 
the three-step nucleon (b) in the case of (ii). In the case of (i), the three-step nucleons are 
formed by the collision of a two-step nucleon with a zero-step nucleon. On the other hand, 
in the case of (ii), the three-step nucleon is produced by a collision of a two-step nucleon 
with a one-step nucleon. In general, the type of collision process in the case of (i), namely 
the collision of a two-step nucleon with a zero-step nucleon, is more frequent in producing a 
three-step nucleon. 

Recently Kawai and his collaborators have developed a semi-classical distorted wave 
model of nucleon inelastic scattering to continuum p0|-p2|. In this model, the first and 
second Born terms are shown to correspond to one-collision and two-collision processes in 
the INC model and are called one-step and two-step processes, respectively. Our above 
definition of one-step and two-step processes is clearly the same as that in this semi-classical 
distorted wave model. We further expect that our iV-step process higher than two-step 
process will be proved to correspond to the iV-th Born term if their semi-classical distorted 
wave model is extended to higher Born terms. 
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III. COMPARISON WITH DATA 



We have studied 12 C(p, p'x) at E p = 200 and 90 MeV, and 58 Ni(p, p'x) at E p = 120 
MeV. The total number of events is 4,000 for 12 C target case while for 58 Ni target case it is 
14,000 for case-1 ctivjv and 3,000 for case-2 and case-3 (Jnn- As stated above, in the case of 
12 C target, we have used only case-1 <Jnn- The data we compare with our calculations are 
taken from Ref. |Q and Ref. |[34] for 12 C and 58 Ni targets, respectively. The range of the 



adopted impact parameter b is < b < 8 fm in the case of 12 C target, while < b < 10 fm 
in the case of 58 Ni target. If the incident proton comes out without making any two-nucleon 
collision, we regard that the event should not be included into inelastic scattering events to 
continuum. In the case of 58 Ni target, we show in Fig. 2 the 6-dependence of the probability 
-Pn.coi.(fr) that the incident proton comes out without making any two-nucleon collision. It is 
to be noted that the total reaction cross section cjr which we discuss later is calculated as 

a R = J q °° 2Mb (l - P n . col .(&)) • (12) 

The calculation of the double differential cross section d 2 a/dQdE p > is made in the fol- 
lowing way; 



dVtdE p i Jo dQdE p r 



,!2Xip ^-dndE p , = ( £ |<p|0 Z! >| 2 W 3 p, (is) 



d{ldE„/ 



= isolated 
protons 



b 



I(P|^)| 2 = (-^) 3/2 exp[--^(p-K J ) 2 



2vr/iV 



2h z u 



where d 3 p = mpdQdE p /, and ( stands for the average value over the events with impact 
parameter b. In this formula, the outcoming protons are expressed by Gaussian wave packets 
with momentum width h,^fv. In the case of protons emitted during the dynamical stage, we 
adopt this width, but in the case of evaporated protons from a fragment of mass number 



Ap, we adopt narrower width by 1/y/Ap, namely hyv/Ap. It is because the center-of-mass 
motion of the fragment of mass number Ap is described by a Gaussian wave packet whose 
momentum width is hy/Apv, from which the standard deviation of momentum per nucleon 
is given by hJv/Ap. In the actual calculation of the angular distribution, we further make 
another modification to | (pl^zj | 2 . It is the cut of the high momentum tail of | (p|0Zi) | 2 - The 
reason of this tail truncation is because the very high momentum part of the distribution 
l(p|^Zi)| 2 has no physical justification. The tail part which we truncate is defined by the 
condition (p — K^) 2 /2% 2 i> > 1.8 2 for dynamical nucleons, and A F (p — Ki) 2 /2h 2 u > 1.8 2 
for evapolated nucleons. The volume of this tail part is about 10 % of the total volume of 
I(P|0 Z! )| 2 . 

In Fig. 3 we give, in the case of the 12 C target, the comparison of the calculated angular 
distribution d 2 a/dfldE p / with the data for various values of outcoming proton energy E p i. 
We see that the reproduction of data is good. It is to be noted that good reproduction 
has been obtained without any adjustable parameter. In the case of E p = 90 MeV, the 
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cross sections at large angles for E p > = 55 and 45 MeV are under estimated. We think that 
the event number 4,000 is still not enough to reproduce large angle cross sections. From 
the study of 58 Ni(p,p'x) shown in Fig. 6, we can guess that we need to increase the event 
number to about 10,000. In Figs. 4 and 5 we show the decomposition of the calculated 
angular distribution into multi-step contributions for the cases of E p = 200 and 90 MeV, 
respectively. We see that for our present region of E p > > 45 MeV, the contribution of 
one-step process is predominant over a wide angular range but yet at large angles we can 
recognize some predominance of two-step process. Especially in the case of highly inelastic 
scattering with E p > = 45 MeV for E p = 200 MeV, the two-step process is dominant over a 
fairly wide angular range and is indispensable for the reproduction of data. 

The comparison of the calculated angular distribution with the data in the case of the 
58 Ni target is given in Fig. 6. We see that the reproduction of data is again good in view of 
the absence of any adjustable parameter in our approach. What is to be noted here is the 
result that almost the same good reproduction of data is obtained for all three choices of 
in-medium cross section a^N, case-1 <Jnn, case-2 <Jnn, and case-3 <Jnn- To state in more 
detail, case-1 <7nn gives slightly better data-fitting for E p > = 100 MeV, while case-3 cr^w 
gives somewhat better fitting for E p > = 40 MeV. 

The decomposition into multi-step contributions of the calculated angular distribution 
is given in Fig. 7. The decomposition of the calculated angular distribution by case-2 cr NN 
is similar to that by case-3 <Jnn- We see that two-step contribution is dominant at some 
large angle region and is indispensable for the reproduction of data. What is to be noted 
is the forward cross section in the case of highly inelastic scattering with E p > = 40 MeV. 
In this case the two-step contribution is much larger than the one-step contribution in the 
forward cross section in all the three cases of adopted in-medium <tnn- It is because the 
quasi-free peak of the one-step process is located in fairly large angle region in the case 
of highly inelastic scattering. Another point to be noted is the contribution of three-step 
process in large angle region. Fig. 7 shows that in both cases of <tnn the data points above 
9 = 100° for Ep> = 100 MeV are reproduced mainly by three-step contributions. Also for 
E p i = 60 MeV, the data points above 9 = 120° have large contributions from the three-step 
process. We clearly see that the two-step contribution is much more important in this case 
of 58 Ni target than in the case of 12 C target, which is of course quite natural. 

In Fig. 7, we indicated with arrows (a) the angle of the quasi-free peak which is given by 



cos -1 J E p i / E p . The peak angles of the calculated one-step angular distributions are however 
seen to be slightly shifted to forward direction. One of the reasons of this shift is the mean 
field effect. If the collision takes place at the point of the potential depth V, the angle of the 



quasi- free peak is given by cos -1 J (E p i — V)/(E P — V). When we take V = —50 MeV, the 
shifts to forward direction of the angles of the quasi-free peak at E p i = 100, 60, and 40 MeV 
are 4°, 9°, and 12°, respectively. We show in Fig. 7 by arrows these shifted angles of 
the quasi-free peak obtained with V = —50 MeV. We see that these shifted angles are now 
much closer to the peak angles of the calculated one-step angular distributions. Another 
possible reason of the angle shift is the refraction effect in the surface region of the target 
which makes the path of one-step nucleon bend to forward direction. 

Recently Shinohara et al. [22] studied this reaction, 58 Ni(p, p'x) at E p = 120 MeV, by 
the use of a semi-classical distorted wave model of nucleon inelastic scattering to continuum 
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2) . Our angular distribution due to the one-step process is very similar to that of Ref . 
|22fl , while the two-step contribution of Ref. [22| is fairly larger than that of our calculation. 
Especially the bigger contribution of the two-step process than the one-step process seen in 
the case of E p i = 60 MeV in Ref. p2l is in disagreement with our result. 



IV. DEPENDENCE ON IN-MEDIUM N-N CROSS SECTION 

We have seen in the previous section that the calculated angular distribution d 2 a/dfldE p/ 
of 58 Ni(p, p'x) at E p = 120 MeV for E p/ > 40 MeV is insensitive to the choice of different 
in-medium cross section aNN- We investigate here the reason of it. 

In Table II we show the calculated reaction cross section a R (Eq. ([12]) ) by the use of 
case-1 cttvat, case-2 &nn, and case-3 (Tnn- We see clear difference between the values of or's 
obtained by different a NN . Reflecting the relation (a A r Ar ) casc _i < (<7 NN ) casc _ 2 < (cr NN ) C£tse _ 3 , 
the resulting a R satisfy (cr K ) case _i < (a R ) casc _ 2 < (o R ) casc _ 3 . As is reflected in Eq. (f[2j), 
the magnitude of a R is determined only by the magnitude of the cross section of the first 
collision of the incident proton. Namely the magnitude of a R does not reflect any information 
of the later reaction process after the first collision of the incident proton. On the other 
hand, the angular distribution for a fixed value of E p i reflects not only first collision of the 
incident proton but also the reaction process after the first collision. Let us compare the 
reaction process due to (o"ATjv) C asc-i with that due to (cr7Vjv) case _3. Since (cr at at ) ca se-3 is larger 
than (o"ArAr) case -i, the probability of the first collision by (crNN)c&se-3 is higher than that by 
(o"AfAf)case-i- But at the same time, the probability of the second and further collisions by 
(c r AfAf)case-3 is also higher than those by (o"jV7v) case _i. Therefore the energy of the outcoming 
protons are much more damped in the average for (o"AT7v) case _3 than for (cx7VAr) case _i. Thus the 
cross section da/dE p i is not necessarily larger for (<JNN)case-3 than for (o"AT7v) case _i, at the E p i 
value where the one-step process makes large contributions. This compensation mechanism 
between the first- collision probability and the energy damping due to multiple collisions 
explains why the calculated angular distribution at fixed value of E p i is almost independent 
of the in-medium N-N cross section within the present three kinds of <Jnn- 

In Table II we have also shown the cross sections o"i_ step of outcoming protons due to the 
one-step process, for three cases of aNN- We see that cri_ stop is smaller for larger aNN- This 
result supports our above argument. If aNN is large, the nucleon after the first collision 
will make further collisions rather than escaping from the target without making further 
collisions. If this effect is bigger than the effect of the large cross section of the first collision, 
tf"i-step becomes smaller for larger aNN- Table II shows that it is actually the case. 

The cross section o"d yn in Table II is the cross section of dynamical protons, namely 
protons emitted during dynamical stage. We see that crd yn is not so much dependent on 
the magnitude of a nn- It is because the cross sections of protons due to higher multi-step 
processes tend to become larger for larger aNN, although o"i_ step is smaller for larger a nn- 

In Table II we also give the calculated values of the cross sections cr evap of evaporated 
protons. Quite naturally cr evap is larger for larger aNN- The cross section a p > in Table II is the 
total cross section of all the outcoming protons, namely oy = a cvap + o"d yn - The dependence 
of a p i on onn is due to that of oWap- 



In the approach with the semi-classical distorted wave model of Refs. [p2| , p0| -|32]| , the 
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imaginary part W opt of the optical potential and the in-medium cross section aNN are im- 
portant ingredients of the model. These W opt and ct^n are mutually intimately related and 
hence the choice of these quantities should be made consistently. It means that the change 
of the magnitude of aNN needs to be made together with the consistent change of W opt . On 
the other hand, in the case of our AMD approach, we need not to introduce W op t and the 
role of W opt is described automatically in the many-body dynamics which includes the two- 
nucleon collision process. Therefore when we change the aNN value in the AMD approach, it 
implies that such effects that are described by W opt are changed consistently in an automatic 
way. This merit of the AMD approach is very advantageous in our present discussion of the 
dependence of d 2 a / dfldE p / on cr NN . It is, however, to be noted that in principle there exists 
mutual relation between a^N and the effective nuclear force. Hence the change of ctjsin is to 
be correlated with the change of the effective nuclear force. In our present study we study 
the dependence of d 2 a j d£ldE p i on a NN by assuming that the adoption of the Gogny force 
as the effective nuclear force is more reliable than the choice of aNN- 

We have seen that the reaction cross section a R depends rather strongly on <Jnn- There- 
fore the comparison of the calculated a R with data is expected to give us important infor- 
mation on a n n- We here should recall that a R is determined by the cross section of the very 
first collision of the incident proton. Hence aNN which we expect to extract from a R data 
is that in the region of the incident energy E p . To the knowledge of the present authors, 
the a R data of the 58 Ni + p system at E p = 120 MeV are not available. However, as we 
discuss below, we have several reasons to expect that cx R ( 58 Ni, 120 MeV) lies in the region 
of 650 ~ 800 mb. In Ref. [35] in which a R of 40 Ca + p, 90 Zr + p, and 208 Pb + p systems are 
experimentally studied in the wide energy range up to 200 MeV, we see that the approximate 
target-mass-number (A T ) dependence and the approximate energy dependence of a R in all 
these three systems are a R oc A 2 ^ 3 and da R /d In (E p /MeV) ~ —170 mb, respectively, in the 
region of 50 MeV < E p < 150 MeV. When we apply this empirical rule of At dependence 
to the data a R ( 40 Ca, 120 MeV) and a R ( 90 Zr, 120 MeV) given in Fig. 12 of Ref. g|, we get 
a R ( 58 m, 120 MeV) « 750 mb. According to Ref. observed a fi ( 58 Ni) at E p = 60.8 MeV 
is 807 mb. When we apply the above empirical rule of energy dependence to this, we obtain 
a fi ( 58 Ni, 120 MeV) « 690 mb. Furthermore in Ref. [0 it is reported that a R ( 59 Co) = 780 
mb at E p = 98.5 MeV and a R (Cu) = 751 mb at E p = 113 MeV. All the above data and 
estimations support that a^( 58 Ni) at E p = 120 MeV is in the region of 650 ~ 800 mb. When 
we accept this range of value for OR( 58 Ni) at E p = 120 MeV, we can extract the conclusion 
that case-1 aNN is most plausible in the region of E p = 120 MeV among the three cases of 
07^, although case-1 a^N gives still rather large a R . 

In Fig. 6 we see that in the case of E p i = 40 MeV the calculated value of d 2 a/dVLdE p i 
is smallest in the case of the case-3 aNN than case-1 and case-2 aNN- This makes the 
reproduction of the observed angular distribution by the case-3 aNN slightly better than 
that by other aNN- Since the contribution of the one-step process is largest in the angle 
region around quasi-free peak, we can understand the situation as follows: In order for the 
proton with E p i Pd 40 MeV after the first collision to escape out of the target, this proton 
should not make second collision. If aNN is large, the probability to make the second collision 
is large and hence the value of d 2 ajd£tdE p i due to one-step process becomes smaller. This 
argument implies that the in- medium cross section a n n in the low energy region around 40 
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MeV is better represented by the case-3 <jnn than other case-1 and case-2 <Jnn- Namely, 
ctjvtv m the low energy region around 40 MeV is suggested to be closer to the free cross 
section to than case-1 and case-2 cr NN . But of course much more investigations are needed 
to have conclusive justification about this point. 

V. SUMMARY 

Proton inelastic scatterings to continuum 12 C(p, p'x) at E p = 200 MeV and 90 MeV and 
58 Ni(p, p'x) at E p = 120 MeV have been studied with AMD (antisymmetrized molecular 
dynamics) which has been successfully utilized for the study of heavy-ion collisions. An- 
gular distributions d 2 a / dVLdE p i for various E p i energies have been shown to be reproduced 
well without any adjustable parameter. It shows the reliability and usefulness of AMD in 
describing light-ion reactions. 

Decomposition of the calculated angular distributions into multi-step contributions has 
been made and two-step contributions have been found to be dominant in some large angle 
region and to be indispensable for the reproduction of data. In the case of highly inelastic 
scattering with E p , = 40 MeV of 58 Ni(p, p'x) at E p = 120 MeV, two-step contributions 
have been found to overwhelm one-step contributions in forward cross sections. Further- 
more in the 58 Ni(p, p'x) case, even three-step processes have been found to make dominant 
contributions in the reproduction of data in the large angle region 6 ^ 120° for E p i = 100 
MeV. 

We have studied the dependence of the calculated angular distributions d 2 a /dVLdE p > for 
various E p i energies on the in-medium N-N cross section <jnn i n the case of 58 Ni(p, p'x) 
at E p = 120 MeV, and have found that the dependence is very weak. The reason of this 
insensitivity of d 2 a / dVLdE p i to ctnn has been clarified by studying also the a^Ar-dependence 
of various kinds of cross sections such as the reaction cross section <tr, the one-step-proton 
cross section o"i_ s t ep , and the evaporated-proton cross section owp. Unlike d 2 a j 'dVtdE p i , these 
quantities, <Jr, 0"i_ stcp , and <T cvap , have been shown to be sensitive to o~nn- For the choice of 
larger (Jnn, &r and o" cvap become larger, while 0"i_ step becomes smaller for larger <tnn- 
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FIGURES 




FIG. 1. Illustrative examples of multi-step process. Nucleon (a) is the incident proton and 

collision points are labeled by collision indices 1, 2, . . . instead of No.l, No. 2, In the figure (i), 

the outcoming nucleons (b) and (c) carry the same set of collision indices {1,2,3} and hence are 
three-step nucleons. In the figure (ii), the exit nucleon (b) is a three-step nucleon since it carries a 
set {1, 2, 3} composed of three collision indices, while the exit nucleon (c) carrying a set {1, 2, 3, 4} 
is a four-step nucleon. 
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FIG. 3. Comparison of calculated angular distribution d 2 a / d£ldE p i of 12 C(p, p'x) with data 
for various E p i energies. Comparisons are made in two cases of incident proton energies E p = 200 
MeV and 90 MeV. 
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FIG. 4. Decomposition of the calculated angular distributions into multi-step contributions 
in the case of 12 C(p, p'x) at E p = 200 MeV. 
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FIG. 5. Decomposition of the calculated angular distributions into multi-step contributions 
in the case of 12 C(p, p'x) at E p = 90 MeV. 
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Decomposition of the calculated angular distributions of 58 Ni(p, p'x) at E p 



120 



MeV into multi-step contributions. The angles indicated with arrows (a) and {(5) are quasi-free 
scattering angles given by cos" 1 E p > / E p and cos" 1 (E p / — V)/ (E p — V) with V = —50 MeV, 
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TABLES 



TABLE I. Values of parameters Ai(E), B\(E), and B2(E) in mb/sr as functions of the incident 
particle laboratory energy E, used in the parametrization of the angular distribution of (Tp n of the 
case-3 <tnn which is taken from Ref. |2{|. See Eq. (|To|). 
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A 1 (E) 


B 1 (E) 


B 2 (E) 


40 
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7.0 
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8.1 
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120 
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240 
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3.6 
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280 


1.8 
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6.0 



TABLE II. Dependence of various kinds of cross sections on the different choices of in-medium 
N-N cross section a^N, case-1 ~ case-3. The notations <tr, 0"i_ s t ep , 0"dyn> oy> an d cr evap stand for 
the reaction cross section, the one-step-proton cross section, the dynamical-proton cross section, 
the total-proton cross section, and the evaporated-proton cross section, respectively. Units are in 
mb. 
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c ) cr dyn 
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d ) oy 


1749 


2259 
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6 ) Ccvap 


722 
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